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We report our recent progress on the nuclear symmetry energy probe, which is called the isobaric yield ratio 
difference (IBD), and its application in neutron density determination in experiments. The results obtained by 
the IBD, from which the isobaric yields in the measured 140A MeV 40.48Ca + Be and **Ni + °Be reactions, 
and the calculated 80A MeV **°*Ca + "C reactions by using a modified statistical abrasion-ablation model, 
show the sensitivity of the IBD to the density differences between reactions. 
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I. INTRODUCTION 


The symmetry energy of asymmetric nuclear matter, 
especially neutron-rich matter, has been an important topic 
in both theory and experiment with the development of 
radioactive nuclear beam (RNB) facilities. The new or up- 
dated RNB facilities can provide nuclear beams with larger 
neutron-richness, which are supposed to promote the study of 
neutron-rich nuclear matter and extend the chart of nuclides 
to the more neutron-rich side. In heavy-ion collisions above 
the intermediate energy, nuclear matter from sub-saturation 
to supra-saturation densities can be produced. The symmetry 
energy of nuclear matter depends on density and temperature. 
In the past decade, though many probes have been proposed, 
there are still difficulties in measuring the symmetry energy 
of nuclear matter at different nuclear densities. 

The isobaric yield ratio (TYR) method was first proposed 
by Huang et al. to extract the symmetry-energy coefficient 
for fragments with small isospins, which also have some tem- 
perature since they are measured in heavy-ion collision [1, 2]. 
In the TYR, some terms deciding the fragment yield cancel 
out, which makes it useful to study the symmetry energy term 
of the fragment with a low temperature. The TYR method 
is extended to study the symmetry energy of fragment with 
larger neutron richness [3—6]. The results obtained by the 
TYR method are compared to a similar method, i.e., isoscal- 
ing, which extracts the symmetry energy of the colliding 
source, showing a large difference between them [7—10]. We 
studied the results of the TYR and isoscaling, and found that 
though the TYR and isoscaling methods both aim at the ex- 
traction of nuclear symmetry energy, the obtained results are 
for different nuclear matter. A new probe was constructed via 
TYR between two reactions, which is called the isobaric yield 
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ratio difference (IBD), and the results of the IBD and isoscal- 
ing were found to be similar because they are both the probes 
for the chemical potential of neutrons and protons in the re- 
actions [11, 12]. Moreover, the IBD probe is also sensitive 
to the density differences between the reactions [13, 14]. In 
this article, we will summarize the recent progress on the IBD 
method. 


Il. THEORETICAL DESCRIPTIONS 


In thermodynamics theory, an equilibrium state is assumed 
for the reaction, which introduces the concept of tempera- 
ture. The temperature affects the yield of fragment. With 
high temperature, light particles like neutrons and protons 
can escape from the binding of hot fragment. The formation 
time of the final fragments, which correspond to the mea- 
sured ones in experiments, is called the “chemical freeze-out” 
(when the fragment ceases to emit light particles, like neu- 
trons, protons or a, and the neutron and proton numbers do 
not change anymore). In general, the yield of fragment in 
heavy-ion collisions above the Fermi energy obeys an expo- 
nential law. For example, in the grand-canonical ensembles 
theory within the grand canonical limit, the yield of a frag- 
ment with mass numbers A and neutron-excess J (I = N-Z 
is the difference between the numbers of neutrons and pro- 
tons) is written as [15, 16] 


o(I, A) = cA" exp{[F(I, A) + Nun + Zitpl/T}, 0) 


where T is the temperature of the system. F(I, A) is the 
free energy of the fragment which depends on T. jun(ftp) is 
the chemical potential of the neutrons (protons) in the equi- 
librium system, which depends on the density and T of the 
systems. In the isoscaling method [17, 18], the difference 
between the un(up) of two reactions can be extracted from 
the isotopic (isotonic) ratio 
o(I , A)2 


Kae TA (2) 


x exp{[N (unz — Hn) + Z(up2 = Hp1)]/T}, 
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with indices 1 and 2 denoting the reactions. In Eq. (2), the 
T of the two reactions is assumed to be the same. For the 
chemical freeze-out fragments this can be true since they both 
have a very long evolution time. For the isotopic ratio, 


ln RU, A)21 = [N (ine 


Uni) 4 Z (up2 Mpi)|/T, (3) 


Z does not change. The isotopic ratio shows scaling phenom- 
ena obeying the linear function. The fitted slope is defined as 
a = (Un2 — Uni)/T. Similarly, for the isotonic ratio, N is 
a constant. The isotonic ratio also shows scaling phenome- 
na obeying the linear function. The fitted slope is defined as 
B = (Up2 — Hpı)/T. a and £ are used to extract the sym- 
metry energy of the system. But the symmetry energy in the 
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isoscaling method is an indirect result obtained from the yield 
of fragments. 

Starting from Eq. (1), assuming the isobars differing 2 units 
in J have the same temperature, the isobaric yield ratio can be 
defined as 


o(I + 2, A) 
a(I, A) (4) 
= exp{[F (I + 2, A) z HL, A) + bn — Mp|/T}. 


R(I4+2,1,A)o1 = 


To make cancellation of F, we adopt a reference reaction 
system with similar measurements as the isoscaling method 
does, and assume that the free energies of one fragment in the 
two reactions are the same. Defined as the difference between 
the IYRs of two reactions, the IBD is [11, 12] 


Alh R(I + 2, I, A)ai = In Ro(I + 2, I,A) = In Ri (I + 2,1, A) 
= [(Hn2 — Uni) — (Mp2 — Mpi))/T (5) 
= (Apn2i = Apy21)/T = Apz1/T. 


The indices also denote the reaction systems. From the 
definition of the isoscaling parameters, the following correla- 
tion can be obtained 


Al R(I +2,1, A) =a -— B. (6) 


Eq. (5) can also be written as 


Aln R(T + 2,1, A)21 = [(Hn2 — Mp1) — (na — Mpi)|/T 
= (Aunp2 = Alnpı )/T. 
(7) 
The chemical potentials of neutrons and protons are related 
to the densities of neutrons (pn) and protons (pp) respective- 
ly [14, 19] 


Un/T =Inpn, and up/T = ln pp, (8) 


which result in a = ln pn2 — ln pni = In(pn2/pni), and 
B = ln pp2 — ln ppı = ln(pp2/Pp1). We can also write the 
relationship between the IBD results and the densities of the 
reaction systems as [14] 


Aln R(I =F 2, I, Ajai = In(Pn2/Pni) = In(Pp2/Pp1); (9) 


Aln RU + 2,2, A)21 = In(pn2/Pp2) — M(Pn1/Pp1). (10) 


Eqs. (9) and (10) are written according to Eqs. (5) and (7), 
respectively. The densities of protons can be determined by 
experimental measurement. Thus the relationships in Eqs. 
(9) and (10) can be used to determine the neutron density of 
the neutron-rich system. For example, if the reference re- 
action system is a neutron-proton symmetric one, in which 
Pni and ppı can be assumed to be the same and pp2 is mea- 
surable, pn2 can be determined from the IBD results. In 


Eq. (9), if pp2 and pp; can be assumed to be the same (for 
isotopic projectiles) and pp; © pni; AR( + 2,7, A) = 
In(Pn2/Ppi). Similarly, in Eq. (10), the reference reaction is 
neutron-proton symmetric and pp; % Pni, by knowing pp2, 
AR(I+2,I, A)21 = In(Pn2/Pp2) and pn2 can be determined 
from the IBD. 


Ill. RESULTS AND DISCUSSION 


We will show the IBD results obtained from the measured 
140A MeV *°8Ca + °Be and *8°Ni + °Be projectile frag- 
mentation reactions [20]. These reactions were measured by 
M. Mocko et al. at the National Superconducting Cyclotron 
Laboratory (NSCL) at Michigan State University. The mea- 
sured fragments range from Z = 5 to Zpro;, which provides 
high quality data for the test of the IBD methods. The IBD 
results has been previously reported [11]. 

In Fig. 1, the results of the IBD and isoscaling methods 
are re-plotted for comparison [11]. Four groups of reactions 
are analyzed, the **Ca/*°Ca + °Be and “Ni/*®Ni + °Be reac- 
tions of the isotopic projectiles, the °°Ni/*°Ca + °Be reactions 
with the projectiles both being neutron-proton symmetric, and 
the *8Ca/*Ni + °Be reactions with the projectiles both being 
neutron-rich. The reactions using the relative neutron-proton 
symmetric projectile was set as the reference reaction denoted 
by index 1. Theoretically, the relationship between the IBD 
and isoscaling is described in Eq. (6). The fragment belongs 
to one isotopic chain and one isotonic chain, thus it can be 
related to a and 8 according to the respective isotopic and 
isotonic ratio fitting. For the isobaric chains with J from —1 
to 2, the IBD and isoscaling results were found to be similar 
except for some fragments with a large A. The IBD results 
were the same in the J = 0 and 1 isobaric chains for the 
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Fig. 1. (Color online) The comparison between the IBD and isoscaling results (a re-plotted figure from Ref. [11]). The circles and squares 
denote the IBD and isoscaling results, respectively. The lines are just for guiding the eyes to the plateaus in the IBD and isoscaling results. 


48Ca/*Ca reactions. The same phenomena can be found in 
the 7 = 1 and 2 isobaric chains for the “Ni/>8Ni reactions, 
and the 7 = —1 isobaric chains for the °*Ni/*°Ca reactions. 
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Fig. 2. (Color online) The asymmetry of the systems for obtaining 
the IBD results. The labels on the circles denote the reactions used. 


It was also found that the height of the IBD-plateaus in the 
four groups were different, about 2, 1.4, 0.5 and 0.4, respec- 
tively [11]. And the lengths of the plateaus are also different. 
In the IBD of neutron-rich isobars, the plateau can disap- 
pear [11-14]. The IBD results denote the difference between 
the neutrons and protons for the two reactions according to 
Eqs. (9) and (10). If the IBD results change very small, it 
reflects that the difference between the densities of neutrons 
and protons for the reactions are small. The increasing IBD 
with A of the fragment indicates an enlarged difference be- 
tween the densities of the neutrons and protons [21, 22]. Tak- 
ing the IBD results for the **Ca/*°Ca reactions as an example, 


Pp1 and pnı for 40Ca can be assumed to be the same, and as- 
suming that Ppl © Pp2, AR(I + 2, I, A)21 = In(pn2/Pp1)- 
Pn2 Of the 48Ca reaction can be determined from the IBD 
results. For the “8Ca/*Ca reaction, AR(I + 2,1, A)21 = 
In(&*) = 2 which results in a = 7.4. Itis very large 
compared to the isoscaling method. It should be noted that 
48Ca has a much larger neutron excess than *°Ca, and a ~ 
8, of which 8 cannot be omitted in Eq. (6). This makes 
Aln R(I + 2,2, A)21 = 2a , which results in the overes- 
timation of ne Similar explanations can be found for the 
other three reaction groups. It should be noted that the densi- 
ties determined by the IBD results are for the time of freeze- 
out, assuming thermal equilibrium. The relationship between 
the density of the projectile and the system at the freeze-out 
should be further studied in transport models. 


The IBD-plateau was found to be sensitive to differences 
between the densities of the reaction systems. The 
asymmetry of these systems (N/Z)21 are calculated 
and the results are plotted in Fig. 2. The asym- 
metry of the reaction systems (N/Z)2, is defined as 
[Not MW Zp+Z)lo/ Np +N Z+]. From high to low, 
(N/Z)o1 is in the order of *8Ca/*Ca, *8Ca/*8Ni, “Ni/*°Ca, 
4NG/>8Ni, *Ca/Ni, and >8Ni/*’Ca. The (N/Z)21 of the 
48Ca/“Ni and ~8Ni/*°Ca reactions are very similar. In 
Fig. 3(a), the IBD-results for the J = 0 isobaric chain are 
re-plotted. From high to low, the IBD-plateau is in the order 
of *8Ca/Ca, Ca/>®Ni, “Ni/*Ca, “Ni/*8Ni, *8Ni/*Ca, and 
48Ca/“4Ni. The IBD-plateaus for the *8Ni/*°Ca and *8Ca/“Ni 
reactions are very similar. The order of the (N/Z)2: and 
the IBD-plateau are consistent, except for the *8Ni/#°Ca and 
48Ca/4Ni reactions. The same phenomena can also found in 
the IBD results for the J = 1 isobaric chain in the reactions 
as plotted in Fig. 3(b). 

To test the sensitivity of the IBD-plateau to the density 
difference between the reactions, we calculate the fragment 
yield in the 80A MeV ?8-5?Ca + !°C reactions (only the even- 
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Fig. 3. (Color online) The IBD results for the J = 0 (a) and J = 1 
(b) isobaric chain. 


A calcium isotopes are used) by using a modified statistical 
abrasion-ablation (SAA) model [6, 13, 14, 21-27]. The sta- 
tistical abrasion-ablation model will not be described in this 
work. Readers can refer to Refs. [21-27]. The fragments pro- 
duced in the Ca + !C reaction is set as reaction 1, which is 
the reference reaction. The IBD results for the isobaric chains 
from I = —1 to 2 are plotted in Fig. 4. Since the analyzed 
fragments are the prefragments which are determined by the 
neutron and proton densities, as well as the nucleus-nucleus 
reaction cross sections, the yields of prefragments are sen- 
sitive the density changes of the projectile [13]. It is found 
that, as the projectile (*Ca, with X denotes the mass num- 
ber of projectile) becomes more neutron-rich, the IBD results 
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Fig. 4. (Color online) The IBD results for the prefragments calculat- 
ed by using a modified statistical abrasion-ablation model. 
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Fig. 5. (Color online) The averaged value of IBD-plateaus for the 
isobaric chains from J = —1 to 2. The inserted figure show the 
difference between the neutron density distributions of “Ca and “Ca 
with A In pn = In pn (“Ca)— In pp (*°Ca). 


increases regularly. The IBD obtained from the prefragments 
are sensitive to the difference between the neutron density of 
the * Ca and *°Ca reactions if we assume that the proton den- 
sities in * Ca are the same or vary very little [13, 14]. 

We calculated the average values (< A/T >) of the IBD- 
plateau for the isobars plotted in Fig. 4. The < Ay/T > are 
plotted as a function of J in Fig. 5. It is shown that for the 
relative symmetry fragments, < Ay/T > changes very small 
with J. In the inserted figure of Fig. 5, the difference between 
the neutron densities (A ln pn) of *Ca and *°Ca is plotted. 
When r <~ 2 fm, Aln pn are flat and change very small. 
It has been illustrated that the yield of fragments with small 
mass numbers are mostly produced in the central collisions 
[21, 22, 26, 27]. The flat distribution of A In pn between the 
core part of the calcium isotopes can account for the phenom- 
ena of the IBD-plateau. The enlarged A In pn with r (r >~ 2 
fm for different calcium isotopes) can also account for the in- 
creasing of Ay/T with A, as well as the disappearance of the 
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IBD-plateau in the neutron-rich isobars [11-14]. 


From both the experiment results and the calculated results 
by the statistical abrasion-ablation model [11-14], itis shown 
that the IBD probe is sensitive to the density difference be- 
tween reactions. The IBD probe provides a new method to 
determine the neutron density in heavy-ion collisions. But 
as we have noted, the density is for the time of the chemi- 
cal freeze-out, assuming there is a thermal equilibrium in the 
reaction. Actually, the process of the reaction is dynamic 
and changes with reaction time. The IBD results should be 
verified by a dynamic evolution process. Previously we have 
reported the IBD result for the **Ni + !?C reaction calculated 
by the anti-symmetric molecular dynamical (AMD) model to 
verify the sensitivity of the IBD to density [11]. In dynamical 
model calculations, the evolution of the IBD results with time 
can be studied. We have finished the calculations and the IBD 
analysis of the 5864Ni + !?C reactions and hope to report the 
AMD results in near future to see the IBD for primary frag- 
ments and fragments after the decay process. Moreover, the 
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